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EFFECTS  OF  CURVATURE  AND  DILUTION  ON  UNSTEADY,  PREMIXED, 
LAMINAR  FLAME  PROPAGATION 

Introduction 

In  this  paper,  we  present  and  disouss  time -dependent  calculations  of 
one-dimensional  laminar  flames  in  stoichiometric  hydrogen-oxygen  mixtures 
diluted  with  nitrogen.  Ve  use  these  calculations  to  study  the  effects  of 
curvature  (or  stretch)  and  dilution  on  flame  propagation  in  pre-mixed 
gases.  The  results  presented  also  have  important  applications  to  laminar 
flame  quenching  and  flammability  limits  in  the  absence  of  external  heat 
sinks. 

Previous  work  on  propane-air  mixtures  by  Strehlow  (1984)  shows  that 
rich  mixtures  exhibit  a  higher  flame  velocity  at  small  radii  than  at  large 
radii,  while  lean  mixtures  show  a  lower  flame  velocity  at  small  radii  than 
at  large  radii.  In  all  cases,  for  large  enough  radii  the  flame  velocity 
relaxes  to  the  appropriate  planar  flame  velocity.  Strehlow  attributes 
these  effects  to  preferential  diffusion  of  the  lighter  species  and  this 
has  been  verified  by  Frankel  and  Shivashinsky  (1983)  using  an  asymptotic 
a.  lysis.  More  recently,  Law  (1984a;  1984b)  used  an  asymptotic  analysis 
to  explain  the  behavior  of  stretched  flames  in  rich  and  lean  mixtures 
of  methane  and  propane  in  air.  These  theoretical  analyses  are  for  either 
fuel-rich  or  fuel-lean  conditions  and  use  a  simple  phenomenology  to 
represent  the  chemical  kinetics. 

This  paper  presents  calculations  of  planar  and  spherically  expanding 
flames  in  stoichiometric  hydrogen-oxygen  mixtures  diluted  with  nitrogen. 
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The  results  presented  here  ere  obtained  from  numerical  simulations,  which 
use  a  time-dependent,  one-dimensional,  Lagrangian  model  (Kallasanath  et 
al.  1982a).  This  numerical  model  was  developed  specifically  to  study  the 
initiation,  propagation  and  quenching  of  laminar  flames.  In  addition  to 
using  a  detailed  chemical  reaction  mechanism,  the  model  inoludes  the 
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effects  of  molecular  diffusion,  thermal  conduction  and  thermal  diffusion 
of  the  individual  species  considered.  Because  of  the  level  of  detail 
incorporated  in  the  model,  the  spatial  structure  and  temporal  evolution 
of  the  flame  structure  can  be  highly  resolved. 

The  Numerical  Model 

The  numerical  model  solves  the  time-dependent  conservation  equations 
for  mass,  momentum  and  energy  (Williams  1965;  Oran  and  Boris  1981a).  The 
model  has  been  used  for  a  variety  of  flame  studies,  including  calculations 
of  minimum  ignition  energies  (Oran  and  Boris  1981b;  Kallasanath  et  al. 
1982b),  quench  volumes  (Kallasanath  et  al.  1982b)  and  burning  velocities 
(Kallasanath  et  al.  1982c).  The  model  has  a  modular  form  and  permits  a 
wide  variety  of  geometric,  initial,  boundary,  and  time  varying  energy 
input  conditions.  The  algorithms  representing  the  various  chemical  and 
physical  processes  are  Integrated  separately  and  then  asymptotically 
coupled  by  time-step  splitting  techniques  (Oran  and  Boris  1981a).  The 
convective  transport  is  solved  by  the  algorithm  ADINC,  a  Lagrangian  fluid 
dynamic  algorithm  which  solves  implicitly  for  the  pressures  (Boris  1979) • 
The  method  gives  an  accurate  representation  of  material  interfaces  and 
allows  steep  gradients  in  species  and  temperature  to  develop  and  be 
maintained.  In  addition  to  considering  the  thermal  conduction  and 
moleoular  diffusion  processes  in  detail,  the  model  also  includes  thermal 
diffusion.  The  chemioal  interactions  are  described  by  a  set  of  nonlinear, 


coupled  ordinary  differential  equations  which  are  solved  using  a  fully 
vectorized  version  of  the  selected  asymptotic  integration  method  CHEMEQ 
(Young  and  Boris  1977;  Young  1980). 


For  this  study  we  have  used  the  hydrogen-oxygen  reaction  scheme 
(Buries  and  Oran  1981)  which  involves  the  eight  reactive  species  H 2,  02, 

H,  0,  OH,  HgO,  HC  ;,  HgOg  and  diluent,  which  is  chosen  to  be  nitrogen.  The 
thermo chemical  properties  of  the  various  species  Involved  are  taken  from 
the  JANAF  tables  (Stull  and  Prophet  1971).  The  chemical  reaction  mech¬ 
anism  (given  in  Kailasanath  et  al.  1982a  and  1982b)  has  been  extensively 
tested  and  shown  to  give  good  results.  Buries  and  Oran  (1981)  showed  that 
the  results  computed  with  this  mechanism  compare  well  with  experimentally 
observed  induction  times,  second  explosion  limits  and  the  temporal  behav¬ 
ior  of  reactive  species.  Oran  et  al.  (1982)  have  shown  that  the  mechanism 
coupled  to  a  convective  transport  algorithm  gives  good  results  in  the 
simulation  of  the  conditions  behind  a  reflected  shock.  The  laminar 
burning  velocities  calculated  using  the  mechanism  are  in  agreement  with 
experimental  data  (Kailasanath  et  al.  1982c). 

For  the  calculations  presented  below,  the  model  was  configured  with 
an  open  boundary  at  one  end  to  simulate  an  unconfined  system.  Most  of  the 
calculations  were  done  in  a  spherically  symmetric  one-dimensional  geometry 
Other  calculations  were  done  assuming  a  planar  configuration.  All  the 
calculations  are  for  stoichiometric  mixtures  of  hydrogen  and  oxygen  at  an 
initial  temperature  and  pressure  of  298  K  and  1  atmosphere,  respectively. 
The  amount  of  nitrogen  was  varied,  thereby  varying  the  dilution. 

Results  and  Discussion 
Estimation  of  Burning  Velocity 

For  either  thin  or  planar  flames,  the  instantaneous  normal  burning 


velocity  can  be  calculated  froa  the  flame  velocity  if  we  know  the  velocity 
of  the  unburnt  gases  ahead  of  the  flaae.  For  planar  flames,  the  velocity 
of  the  unburnt  gasea  ahead  of  the  flaae  ia  constant,  aa  shown  in  Fig.  1, 
where  the  apatial  variation  of  the  flow  velocity  and  the  temperature 
aero aa  the  flaae  la  ahown  for  a  hydrogen-air  mixture.  Hence  the  burning 
velocity  can  be  unaabiguoualy  determined  aa  the  difference  between  the 
flame  velocity  and  the  flow  velocity  ahead  of  the  flame. 


V.  a  V#  -  V„ 
b  f  flow 


For  a  thin  flaae,  a  similar  definition  ia  adequate.  However,  when  the 
flame  haa  a  finite  thickneaa  and  ia  curved  (aa  for  a  spherically  expanding 
flame),  the  appropriate  definitions  for  the  location  of  the  flame  front 
and  the  fluid  velocity  of  the  unburnt  gaaea  are  ambiguous.  This  can  be 
seen  in  Fig.  2,  in  which  the  spatial  variation  of  the  flow  velocity  and 
the  temperature  across  a  spherical  flaae  in  a  hydrogen-air  mixture  is 
shown.  The  fluid  velocity  reaches  a  maximum  within  the  flame  and  then 
decreases  ahead  of  the  flue. 

For  the  flames  studied  in  this  paper,  two  reference  fluid  velocities 
have  been  chosen  for  estimating  the  burning  velocity.  One  of  them  is  the 


maximum  fluid  velocity  (V  )  in  the  system  and  the  other  is  the  fluid 


velocity  of  the  unburnt  gases  (V^qq)  corresponding  to  the  first  location 


ahead  of  the  flame  with  temperature  of  300  K.  The  lower  estimate  for  the 


burning  velocity  is  obtained  as  (V.  -  V  „  )  and  the  upper  estimate  as  (V_ 

r  Sax  r 


-  ¥jq6).  For  a  planar  flame  in  which  the  fluid  velocity  ahead  of  the  flame 


is  constant  and  the  same  as  the  maximum  fluid  velocity,  the  two  estimates 
for  the  burning  velocity  are  identical. 


Effects  of  Curvature 

In  the  first  set  of  calculations,  a  spherically  expanding  flame  in 
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a  stolohiometric  hydrogen-air  (actually  H2:02:N2  /  2:1:4)  mixture  was 
studied.  In  these  simulations,  energy  was  deposited  linearly  over  a  fixed 
period  of  time  at  the  center  of  a  spherically  symmetric  system.  The 
radius  of  energy  deposition  was  larger  than  the  quench-radius  (Kailasanath 
et  al.  1982b)  and  was  held  constant.  We  then  tracked  the  spatial  location 
of  the  flame  kernel  as  a  function  of  time,  and  used  this  to  calculate  the 
apparent  flame  velocity  as  a  function  of  time.  The  results  of  such  a  cal¬ 
culation  are  shown  in  Fig.  3.  The  flame  velocity  (Vf)  initially  decreases 
with  time  until  it  reaches  a  minimum  value,  and  then  it  increases.  The 
figure  also  shows  the  maximum  fluid  velooity  (Tm  )  of  the  system  and  the 
fluid  velooity  of  the  unburn  gases  (V^qq)  corresponding  to  the  first 
location  ahead  of  the  flame  with  temperature  of  300  K.  With  increasing 
radii,  the  difference  between  the  two  fluid  velocities  decreases.  For 
very  large  radii  (not  shown  in  the  figure),  there  is  an  unambiguous  burning 
velocity  since  the  two  fluid  velocities  approach  the  same  value.  This 
estimated  burning  velocity  approaches  both  the  experimental  value  and  the 
value  determined  from  a  separate  planar  one-dimensional  calculation. 

Effects  of  Dilution 

We  also  performed  a  series  of  calculations  in  which  the  amount  of 
diluent  was  varied.  The  results  of  one  such  calculation,  for  the  H2:02:N2 
/  2:1:7  mixture,  is  shown  in  Fig.  4.  Increasing  the  amount  of  diluent 
does  not  change  the  observed  trends:  the  flame  velooity  and  the  burning 
velocity  decrease  with  Increasing  radii,  attain  a  minimum  value  and  then 
increase  with  increasing  radii.  This  calculation  and  the  results  shown  in 
Fig.  5  for  the  H2:02:N2  ^  2:1:10  oase  show  that  it  takes  longer  for  the 
flame  to  reach  the  minimum  flame  velocity  as  dilution  increases.  For  the 
2:1:10  mixture,  we  expect  the  flame  velooity  to  increase  again  at  larger 
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radius  following  the  trend  seen  in  the  other  two  mixtures.  For  this  mix¬ 
ture  the  burning  velocity  corresponding  to  the  minimum  flame  velocity  is 
.30  -  .35  m/s,  where  as  the  calculated  burning  velocity  for  a  planar  flame 
in  the  same  mixture  is  .85  -  .90  m/s.  Thus  the  burning  velocity  of  a 
spherical  flame  corresponding  to  the  minimum  value  of  the  flame  velocity 

is  smaller  than  that  of  a  planar  flame  in  the  same  mixture.  We  have 

observed  this  trend  in  all  of  the  stoichiometric  mixtures  we  have  studied. 

Some  of  the  more  dilute  mixtures  studied  do  not  support  spherical 
flame  propagation.  One  such  case  is  the  H2:02:N2/  2:1:13  mixture,  in 
which  th«r  flame  velocity  does  not  level  off  but  continues  to  decrease 

until  the  flame  dies.  Figure  6  shows  V.,  7  and  V_nn  as  a  function  of 

r  max  juu 

time  for  this  case.  We  varied  the  radius  of  energy  deposition,  the  amount 
of  energy  input  and  the  mode  of  energy  input.  In  all  cases  the  flame  died 
after  propagating  for  a  short  time. 

When  we  modelled  flame  propagation  in  the  planar  geometry,  we 
observed  signlflciant  burning  velocities  for  a  wider  range  of  mixtures 
than  in  the  spherical  geometry.  For  example,  the  burning  velocity  was 
between  .34  and  .40  m/s  for  the  H2:02:N2  /  2:1:13  mixture.  As  mentioned 
above,  this  mixture  did  not  support  spherical  flame  propagation. 

Effects  of  Stretch 

The  observed  time  histories  of  the  flame  velocity  can  be  explained 
on  the  basis  of  stretch  (due  to  flow  divergence)  and  chemical  kinetics. 
First,  consider  the  initial  deceleration  of  Vf  in  spherical  geometry. 
Initially,  due  to  flow  divergence  the  energy  released  in  chemical  reac¬ 
tions  does  not  balance  the  energy  conducted  and  diffused  into  the  unre¬ 
acted  mixture.  Because  of  this,  the  flame  velocity  and  temperature 


decrease  as  the  flame  expands.  The  decrease  in  the  flane  temperature 
decreases  the  energy  release  rate,  which  leads  to  a  lower  flame  velocity, 
and  so  on.  However,  this  prooess  does  not  continue  indefinitely  because 
the  stretch  effect  decreases  with  increasing  radii.  Because  of  this 
effect  a  minimum  flame  velocity  is  observed  in  spherical  flame  propaga¬ 
tion.  As  the  radius  of  the  flame  increases  further,  the  energy  i ^leased 
in  chemical  reactions  is  larger  than  the  energy  conducted  and  diffused 
into  the  unburnt  mixture.  How  the  flame  velocity  increases  with  increas¬ 
ing  radii  until  a  balance  is  attained  between  the  energy  released  by 
chemical  reactions  and  the  energy  conducted  and  diffused  into  the  unburnt 
mixture  ahead  of  the  flame.  Such  a  balance  occurs  for  large  radii  and 
then  the  flame  propagates  as  if  it  were  planar. 

The  above  observations  are  verified  by  the  temperatures  calculated  in 
the  numerical  simulations.  The  time  history  of  the  maximum  temperature  in 
a  spherically  expanding  flame  in  a  hydrogen-air  mixture  is  shown  in  Fig. 

7.  The  maximum  temperature  corresponding  to  the  minimum  velocity  spherical 
flame  is  significantly  lower  than  that  of  a  planar  flame  in  the  same 
mixture.  After  reaching  a  minimum,  the  flame  temperature  increases  and 
tends  towards  the  temperature  corresponding  to  a  planar  flame  'n  the  same 
mixture. 

The  actual  radius  at  which  the  minimum  flame  velocity  occurs  and 
the  radius  at  which  the  effect  of  stretch  becomes  negligible  depend  on  the 
flame  thicknesss.  Since  the  flame  thickness  increases  with  increasing 
dilution,  a  flame  in  a  more  dilute  mixture  will  have  to  propagate  to  a 
larger  radius  in  order  to  overcome  the  effects  of  stretch.  This  explains 
the  observation  made  above  which  showed  that  the  minimum  flame  velocity  is 
reached  at  a  later  time  (and  a  larger  radius)  as  dilution  is  increased. 


Not  only  la  the  aininua  flaae  velocity  reached  later,  but  the  magnitude  of 
the  velocity  and  the  flame  temperature  are  smaller.  Therefore  by  increas¬ 
ing  the  dilution,  one  can  obtain  a  mixture  which  will  be  quenched  by 
endothermic  reactions  before  it  can  overcome  the  effects  of  stretch.  This 
appears  to  be  the  case  with  the  H2:02:N2  /  2:1:13  mixture. 

Summary  and  Conclusions 

We  have  studied  flames  in  stoichiometric  hydrogen-oxygen  mixtures 
diluted  by  nitrogen  in  spherical  and  planar  geometries  using  time-depen- 
dent,  one-dimensional  numerical  simulations.  The  model  includes  a  full 
description  of  the  eight  reactive  species  involved  in  hydrogen-oxygen 
kinetics  plus  detailed  models  for  the  chemical  kinetics,  thermal  conduc¬ 
tion,  and  thermal  and  molecular  diffusion  processes.  These  studies  have 
led  to  a  number  of  Interesting  observations  and  important  questions  on 
flame  behavior. 

We  have  seen  that,  in  a  stoichiometric  hydrogen-air  mixture,  the 
flame  velocity  of  spherically  expanding  flames  first  decreases  with 
increasing  radii,  and  then  increases.  For  large  radii,  the  burning 
velocity  approaches  the  planar  burning  velocity.  These  same  trends  are 
observed  when  the  amount  of  diluent  is  increased.  However,  with  increase 
in  dilution,  the  minimum  flame  velocity  is  reached  at  a  later  time  and  at 
a  larger  radius.  Comparing  these  spherical  geometry  results  to  another 
set  of  calculations  in  planar  geometry  shows  that  the  minimum  burning 
velocity  of  a  spherical  flame  is  less  than  that  of  a  planar  flame  in  the 
same  mixture. 

In  both  planar  and  spherical  geometries,  the  effect  of  increasing 
the  dilution  is  to  lower  the  burning  velocities.  Since  the  burning 


velocity  is  smaller  in  the  spherical  geometry,  the  flame  can  be  exting¬ 
uished  (or  quenched)  with  less  dilution  in  the  spherical  geometry  than  in 
the  planar  geometry.  The  existence  of  a  quench-radius  and  a  minimum 
ignition  energy  in  the  absence  of  heat  or  radical  loss  to  walls  or  other 
external  sinks  has  already  been  demonstrated  (Kailasanath  et  al.  1982b). 
The  results  presented  here  indicate  that  such  "self-quenching"  also 
depends  on  the  geometry  or  curvature. 

The  time  histories  of  the  flame  velocity  are  controlled  by  stretch 
(due  to  flow  divergence)  and  energy  release  by  chemical  reactions.  The 
maximum  temperature  of  a  spherically  expanding  flame  exhibits  the  same 
trends  as  the  flame  velocity.  The  actual  radius  at  which  the  minimum 
flame  velocity  occurs  and  the  radius  at  which  the  effect  of  stretch 
becomes  less  important  depends  on  the  thickness  of  the  flame.  The  thick¬ 
ness  of  the  flame  increases  with  increasing  dilution.  Therefore  a  more 
dilute  flame  will  have  to  propagate  to  a  larger  radius  in  order  to  over¬ 
come  the  effects  of  stretch.  The  minimum  flame  velocity  and  the  flame 
temperature  decrease  with  increasing  dilution.  Therefore  by  increasing 
the  dilution  one  can  obtain  a  mixture  which  will  be  quenched  by  endo¬ 
thermic  reactions  before  it  can  overcome  the  effects  of  stretch.  This 
appears  to  be  the  case  with  the  1: 13  mixture  which  does 

not  support  spherical  flame  propagation. 

It  has  been  shown  experimentally  that  the  addition  of  inert  gases 
to  confined  fuel-air  mixtures  can  cause  the  flammability  limits  to  become 
narrower  (Zabetakis  1965;  Strehlow  1984).  That  is,  by  adding  sufficient 
amounts  of  an  inert  gas  we  can  obtain  a  mixture  which  does  not  support 
flame  propagation.  Our  calculations  support  this  conclusion  even  though 
we  do  not  Include  any  heat  or  radical  loss  to  the  confining  walls.  The 
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actual  Halt a  obtained  might,  however,  depend  on  parameters  such  as  the 
method  of  Initiation,  the  duration  of  energy  deposition  and  the  geometry 
of  the  system.  By  systematically  varying  these  parameters  as  well  as  the 
stoichiometry,  the  numerical  simulations  can  be  used  to  help  answer 
questions  regarding  the  existence  of  fundamental  flammability  limits. 

Acknowledgments 

This  work  was  sponsored  by  the  Office  of  Naval  Research  through  the 
Naval  Research  Laboratory.  The  authors  would  like  to  thank  Dr.  Jay  Boris 
for  the  many  fruitful  discussions  we  have  had  on  the  ignition,  quenching 
and  propagation  of  flames. 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 

POSITION  (mm) 

Fig.  2  The  flow  velocity  and  temperature  profiles  at  a  particular 
time  in  a  spherically  expanding  flame  in  a  H2:02:1*2/2: 
mixture. 
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Fig.  3  Tine  history  of  the  propagation  of  a  spherically  expanding  flame 
in  a  H2:02:H2/2: 1:4  mixture.  The  flame  velocity  is  denoted  Vf, 
the  maximum  fluid  veloolty  is  VMy,  and  the  velooity  of  the 
first  position  ahead  of  the  flame  with  a  temperature  of  300  K 
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